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Abstract 



Quarkonia have been studied in different collision system and energy in order to un- 
derstand the effects of the hot and dense medium created in heavy-ion collisions. 
CMS is well suited to measure quarkonia decays to muons given the muon identifica- 
tion and charged particle tracking capability. We report here prompt, non-prompt ]/'ip, 
and Y production measured by the CMS experiment in pp collisions at = 7 TeV. 
In addition, the ]/ip and Y production in PbPb at ^/sNN = 2.76 TeV and pp col- 
lisions at the same per nucleon energy are measured and compared. Prompt and 
non-prompt J/i/? contributions are separated for the first time in heavy-ion collisions, 
as is the ground from the excited states in the Y family. Suppression in PbPb at 
^/snn = 2.76 TeV is quantified for prompt ]/tp, B — > ]/ip, and Y(1S), as well as the 
relative suppression of Y(2S+3S) compared to Y(1S). 
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Quarkonia are especially relevant for studying the quark gluon plasma (QGP) since they are 
produced at early times and propagate through the medium mapping its evolution. In par- 
ticular, ]/tp in heavy-ion collisions was suggested to be a promising probe as the deconfined 
medium should screen the two quarks leading to a suppression of its production[lJ. It has 
been studied at different energies and with different collision systems without yet giving a 
fully understood global picture [2-5|. Measuring the charmonium production at the LHC en- 
ergies in PbPb collisions will help constrain predictions, in particular those with a large recom- 
bination probability for prompt ]/tps. Indeed, the abundance of charm quarks in the medium 
could lead to a strong production enhancement at LHC energies [61. In addition to charmo- 
nium precision studies, the LHC center-of-mass energy allows copious Y production in PbPb 
collisions. Detailed measurements of bottomonia will help characterize the dense matter pro- 
duced in heavy-ion collisions complementing the measurements accessible at RHIC energies. 
The full spectroscopy of quarkonium states has been suggested as a possible thermometer for 
the QGP 0. 

This paper first reviews CMS ]/xp and Y cross-section measurements in pp collisions at a/sjvn = 
7 TeV, which allow precision studies of quarkonia production, as well as at a/sjvn = 2.76 TeV, 
which will serve as a reference for the observation of hot nuclear effects in PbPb at the same 
energy. CMS is able to distinguish non-prompt ]/rp from prompt ]/tp in both pp and PbPb 
collisions. The nuclear modification factor (Raa) of prompt, non-prompt ]/ip, and Y(1S) in PbPb 
is measured as a function of transverse momentum {pj), rapidity (y) and number of nucleons 
participating (Npart) in the collision [8J. Finally, the relative suppression of the excited states 
compared to the ground state is quantified [9J. 

A detailed description of the CMS detector can be found in iTO\ . Its central feature is a su- 
perconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within 
the field volume are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter, 
and the brass /scintillator hadron calorimeter. Muons are measured in gas-ionisation detectors 
embedded in the steel return yoke. In addition, CMS has extensive forward calorimetry, in 
particular two steel/ quartz-fiber Cerenkov hadron forward (HP) calorimeters, which cover the 
pseudorapidity range 2.9 < |f/| < 5.2. 

In this paper, quarkonia are identified through their dimuon decay. The silicon pixel and strip 
tracker measures charged-particle trajectories for the range \f]\ < 2.5. The tracker consists of 
66M pixel and lOM strip detector channels, providing a vertex resolution of ~ 15 fim in the 
transverse plane. Muons are detected for the |?/| < 2.4 range, with detection planes based on 
three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. Due to 
the strong magnetic field and the fine granularity of the silicon tracker, the muon transverse 
momentum measurement based on information from the silicon tracker alone has a resolution 
between 1 and 2% for a typical muon in this analysis. CMS is therefore very well suited to 
measure dimuons. In pp, the resolution obtained measuring ]/xp in < 0.5 is 20 MeV/ c^, and 
67 MeV/c^ for the Y(1S) in | < 1. 

J/t/^s can be classified into two types depending on whether they come from the primary vertex 
(prompt ]/tps) or are produced from decays of B mesons (non-prompt ]/rps). Prompt ]/tps group 
direct ]/xp production and ]/ips coming from the feed-down of higher states such as xp' and Xc- 
Non-prompt ]/tps are produced at a distance Lxy from the primary vertex and can therefore be 
separated from the prompt contribution if the resolution of the detector is good enough. This is 
done in CMS by reconstructing the vertices and making a 2-dimensional simultaneous fit 
of the invariant mass distribution and the pseudo-proper decay length, ly^ = Lxy ^ (see [TT] 
for details). 
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CMS has measured prompt and rion-prompt ]/xp cross-sections as a function of pr in different 
rapidity bins [llj. Fig. [T] (left) illustrates the inclusive invariant mass distributions measured 
with Cint = 40 nb"^ 2010 pp data at = 7 TeV over \y\ < 0.5 with the clear ]/xp and the xp' 
peaks. On Fig. [l] (right), the B — > ]/xp CMS cross-section measurement with C^t = 314 nb^^ 
is overlaid with different predictions which reproduces well the data down to low pj in the 
forward region (1.6 < \y\ < 2.4). 
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Figure 1: Left : invariant mass distribution showing the ]/xp and xp' peaks for \y\ < 0.5. Right 
: B — > ]/xp cross-section as a function of pr for 1.6 < \y\ < 2.4 published in [11] compared to 
theoretical predictions. 



CMS is able to identify the three Y states as shown on Fig. |2] (left) and has measured the Y(1S), 
Y(2S) and Y(3S) cross-sections as a function of pr [12] . Fig. |2] (right) shows that the shape 
of the Y(1S) cross-section is well reproduced by PYTHIA but not the normalization which is 
overestimated by a factor two. 
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Figure 2: Left: Y invariant mass distribution for \r]^'\ < 1. Right : Y(1S) cross-section as a 
function of pr for \y\ <2 fT2] compared to Pythia. 



In March 2011, CMS recorded a little more than C^t = 220 nb"^ pp events at ^/s = 2.76 TeV. 
This data is used as a reference for the PbPb measurement. In November 2010, CMS recorded 
-^int = 7.28 jdb^^ of PbPb events, leading to about the same amount of quarkonia statistics as 
the reference pp run at 2.76 TeV. Both data sets have been analyzed following similar condi- 
tions ||8[[9l: (1) events are selected by the CMS two-level trigger keeping any dimuon activity 
in the muon chambers, (2) offline muon reconstruction is seeded with ~ 99% efficient tracks 
in the muon detectors, which are then matched to tracks reconstructed in the silicon tracker by 
means of an algorithm optimized for the heavy-ion environment |[T3l[T4l , (3) the same analysis 
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procedure is followed for the offline selection with very loose criteria. Signal extraction is based 
on the procedures in CMS 7 TeV publications for the signal extraction ||TT1[T2| . 

Using both data sets, the production measured in PbPb collisions is compared to expectations 
from an independent superposition of nucleon-nucleon collisions typically expressed in terms 
of the nuclear modification factor: 



R 



AA 



Cpp Npbpb(QQ) _ £pp 
TaaNmb Npp(QQ) epbPb' 



(1) 



Here Taa is the nuclear overlap functiorj^ Cpp is the pp luminosity, Nmb is the measured num- 
ber of equivalent minimum bias events in PbPb, ^J!"^},^^ is the raw yield ratio, and the 

^ Npp(QQ) ■' '^PbPb 

multiplicity dependent fraction of the efficiency ~ 1-17 for the most central bin). 



Trigger, reconstruction and selection efficiencies of muon pairs are estimated using quarkonia 
Pythia signal embedded in heavy-ion PbPb events generated by Hydjet [18J. These events 
were processed through the trigger emulation and event reconstruction chain. The final effi- 
ciency corrections correspond to the fraction of reconstructed signal passing all the analysis 
selections with respect to the generated signal. Fig. |3] (left) shows this dimuon efficiency as a 
function of Npart for Y (diamonds), prompt (circles) and non-prompt (stars) ]/tp. The individual 
components of the MC efficiency are cross checked using a tag-and-probe technique applied to 
data IITTI . The method consists of fitting the ]/xp candidates, with and without applying the 
probed selection on one of the muons. For example, the dimuon trigger efficiency is probed by 
testing the trigger response from a single-muon triggered sample. Fig. |3] (right) shows that the 
single muon efficiency for the dimuon trigger used in the analysis as a function of measured 
in data (95.1 ± 0.9)% (squares) is compatible with MC (circles). 
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Figure 3: Left : Y (diamonds), prompt (squares) and non-prompt (stars) ]/ip MC efficiency as a 
function of pj in PbPb collisions. Right : single muon efficiency for the double muon trigger 
used as a function on p 7 in MC (circles) and data (squares). 



The CMS detector performs very well in heavy-ion environment such that the good momen- 
tum resolution can be used to separate non-prompt from prompt ]/rp as in pp, making use of 
the distance between the non-prompt vertex and the primary vertex. An example of the 2- 
dimensional fit in PbPb collisions is shown on Fig. |4] for j/tp with pj > 6.5 GeV/ c. For more 
details see ||8l|20l|2Tl. 

^Ratio of the number of binary nucleon-nucleon collisions N„,;; calculated from a Glauber model of the nuclear 
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Figure 4: Pseudo-proper decay length (left) and invariant mass distribution (right) for ]/ip with 
Pt > 6.5 GeV/ c in PbPb collisions.The black dots are the data, the dotted line filled in blue the 
background, the dashed line filled in red the background and non-prompt contribution, and 
the black straight line the total fit. 



For the first time, secondary J/t/^ RA/^ is measured in heavy-ion collisions. Fig. |5] (left) illustrates 
B-meson suppression through their ]/tp decays through the Raa as a function of Npart- Raa = 
0.36 ± 0.08 (stat) ± 0.03 (syst) in the 20% most central collisions. This could be a hint of b-quark 
energy loss. The level of suppression is of the same order of magnitude as charged hadrons as 
observed on Fig. |5] (right) where the non-prompt ]/ip Raa is plotted as a function of pj for 0- 
20% while the bosons and charged hadrons are presented as a function of the transverse mass 
for 0-10% I1T9I. 
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Figure 5: Left: Non-prompt ]/xp Raa in three centrality bin 0-20% and 20-100% in closed sym- 
bols and 0-100% in open symbols. Right: Raa vs. mj for Z (squares), isolated photons (closed 
circles) and charged hadrons (open circles) compared to the secondary ]/ip measurement as a 
function of px (stars). 



Fig. |6] shows the prompt ]/ip Raa (filled squares) as a function of px, y and Npart- A factor 
three suppression is observed for the two pj bin. CMS points are compared to measurements 
at a/snn = 200 GeV from PHENIX [22J at mid- (open squares) and forward (open circles) 

collision geometry I15II16I and the inelastic nucleon-nucleon cross section a^_^^ = (64 ± 5) mb at ys = 2.76 TeV fVT] 
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rapidity for lower pjs, and from STAR ||23| up to px = 8 GeV/ ^ The tendency of high pj 
]/ip's to survive at RHIC is not seen at the LHC. Furthermore, CMS measures less suppression 

at forward rapidity for high pj ]/'ip. One should remember that the x probed with (pj''') = 
10 GeV/c by CMS over \y\ < 2.4 are xi ~ 0.02 and X2 ~ 5 ■ 10^"^ . Therefore, anti-shadowing 
could play a role in the suppression observed and could contribute to seeing an opposite trend 
than PHENIX as a function of y, or an increase of the Raa when going to low pj and more 
forward regions as for ALICE measurements [24J. For pj > 3 GeV/c and 1.6 < y < 2.4, CMS 
measures Raa = 0.39 ± 0.06(stat.) ± 0.03(syst.). Finally, in the 10% most central collisions, CMS 
observes a factor five suppression much greater than measured by STAR. 
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Figure 6: CMS prompt ]/tp Raa measurement (filled squares) compared to PHENIX mid (open 
squares) and forward (open circles) rapidity measurement and STAR higher pj measurement 
(stars) as a function of pj (left), y (center), and Npart (right). 



CMS is able to disentangle the Y(1S) contribution from the higher states in PbPb as in pp 
collisions. Fig. |7| compares the Y invariant mass distribution at = 2.76 TeV in pp (left) 
and PbPb (right) collisions, for > 4 GeV/ c. The higher state contribution relative to the 
ground state is strikingly smaller in PbPb collisions. In order to quantify this suppression, an 
extended unbinned maximum likelihood simultaneous fit to the pp and PbPb mass spectra is 
performed, following the method described in (H, using the parameters detailed in [25 J. The 
ratio of Y(2S + 3S)/Y(1S) in PbPb and pp benefits from an almost complete cancellation of 
possible acceptance and /or efficiency differences among the reconstructed resonances. The 
double ratio obtained is 

Y(2S + 3S)/Y(lS)|pbPb _ +0.19 / . . ^ I n na (^^^. ) (2) 
Y(2S + 3S)/Y(lS)|pp " -0.15 (statj ±0.03 (syst), (2) 

where the systematic uncertainty (9%) arises from varying the lineshape in the simultaneous fit, 
thus taking into account partial cancellations of systematic effects. Finally, using an ensemble 
of one million pseudo-experiments generated with the signal lineshape obtained from the pp 
data. Fig. [7| (left), the background lineshapes from both data sets, and a double ratio (Eq. |2| 
equal to unity within statistical and systematic uncertainties (absence of a suppression), the 
probability of finding the measured value of 0.31 or a downward fluctuation is estimated to be 
0.9%., corresponding to 2.4 sigma in a one-tailed integral of a Gaussian distribution. 

The Y(1S) suppression has a been studied as a function of pj, y and centrality as shown on 
Fig. |8j A suppression by a factor ~ 2.3 is observed for low pj. This seems to disappear for 
Pj > 6.5GeV/c. The rapidity dependence indicates a slightly smaller suppression at forward 



^PHENIX and STAR measurements are inclusive measurements but the contamination from secondary ]/tp is 
expected to be small at RHIC energies. 
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Figure 7: Invariant mass distribution measured in pp (left) and PbPb (right) collisions at 
2.76 TeV for > 4 GeV/c. 



rapidity. In both cases however, the statistical uncertainties are too large for any strong con- 
clusions. In addition, Y(1S) are suppressed by a factor two in 0-10% central collisions. The 
CMS measurement is compared to STAR Y(l+2+3S) preliminary result in AuAu collisions at 
y^SjvN = 200 GeV ||26|| showing a suppression of the same order of magnitude but with large 
uncertainty. 
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Figure 8: Y(1S) Raa as a function of pj (left), y (middle), and Npart (right), compared to STAR 
inclusive preliminary measurement (stars) for the latter. 



In summary, this paper first presented highlights of CMS quarkonia measurements in pp col- 
lisions at a/s = 7 GeV recorded in 2010. The high statistics accumulated at the LHC allows to 
perform details studies that further constrain production mechanisms. CMS has performed the 
measurements of the prompt and non-prompt J/t/;, as well as of the Y(1S) and Y(2S+3S) mesons 
via their decay into pairs in PbPb and pp collisions at 2.76 TeV. Prompt J/i/? has been 

separated from non-prompt ]/ip for the first time in heavy-ion collisions. A strong suppression 
of prompt J/t/J with pj > 6.5 GeV/c is measured in central collisions, and already in periph- 
eral collisions, showing a clear dependence with centrality. Non-prompt J/i/?, though strongly 
suppressed, show no strong centrality dependence within uncertainties. This is the first hint of 
b-quark energy loss in the hot medium. Furthermore, Y(1S) are suppressed by 40% in the 20% 
most central collisions. The comparison of the ratios of Y(nS)-states in pp and PbPb collisions, 
taken at the same center-of-mass energy, is consistent with the partial disappearance of the 
higher states with respect to the ground state in the PbPb collisions. Those two observations 
could indicate that the Y(1S) suppression is due to the melting of the excited states only in PbPb 



7 



collisions. Measuring the amount of suppression caused by shadowing through pA collisions 
together with more precise measurements will be crucial for interpreting what the melting is 
due to. 
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